Dense Regular Packings of Irregular Non-Convex Particles 
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We present a new numerical scheme to study systems of non-convex, irregular, and punctured par- 
ticles in an efficient manner. We employ this method to analyze regular packings of odd-shaped bod- 
ies, not only from a nanoparticle but also both from a computational geometry perspective. Besides 
determining close-packed structures for many shapes, we also discover a new denser configuration 
for Truncated Tetrahedra. Moreover, we consider recently synthesized nanoparticles and colloids, 
where we focus on the excluded volume interactions, to show the applicability of our method in the 
investigation of their crystal structures and phase behavior. Extensions to the presented scheme 
include the incorporation of soft particle-particle interactions, the study of quasicrystalline systems, 
and random packings. 

PACS numbers: 61.66.-f, 61.46.Df, 82.70.Dd, 02.60.-x 



The synthesis of colloids and nanoparticles has ad- 
vanced enormously over the last decade [l,-^. Currently 
it is not only possible to synthesize spherical particles, 
but also a wide variety of convex faceted shapes, such 
as tetrahedra, cubes and octahedra ^ W . Perhaps the 
most remarkable advancement in synthesis techniques is 
the capability to create with high precision and repro- 
ducibility non-convex, irregular, and even punctured par- 
ticles, e.g., colloidal caps [5 , tetrapods [2 , octapods [6], 
and nanostars [3 . Along with the increased availabil- 
ity of complex shapes, there is the concurrent increase 
in the study of their self-assembly into liquid [7 , amor- 
phous [8, 9 , and ordered phases [ilOJ. Especially the or- 
dered (quasi) crystalline structures p!QUT3] of odd-shaped 
particles, as well as their material properties, have re- 
ceived a lot of attention in the past years. Interestingly, 
the research into dense packings of particles is not re- 
stricted to materials science, as it is also connected to 
fields as diverse as discrete geometry and number the- 
ory [14-16^ , computer science and biophysics [17 . 

Predictions obtained from computer simulations on 
the phase behaviour and the self- assembled structures 
of these particles have been essential in guiding experi- 
mental studies and in answering fundamental mathemat- 
ical questions on the packing of particles. Convex ob- 
jects such as spheres [3 HI] and ellipsoids [9], as well as 
(semi) regular [TsHlSl [TSHST] and space-filling [27 solids 
have been the subject of intense ongoing investigation. 
However, ordered structures comprised of irregular non- 
convex particles have hardly been studied by simulation 
because of the numerical challenges in implementing ex- 
cluded volume interactions for such systems, due to the 
complex particle shape and the additional rotational de- 
grees of freedom. Only recently were the first attempts 
made to study such systems, namely for superdisks and 
superballs [23]. 



In this Letter, we present a novel numerical method 
to study systems of non-convex irregular poly topes, col- 
loids, and nanoparticles. We employ this method to es- 
tablish a rigorous lower bound ^ib to the volume fraction 
of their densest regular packing and to predict candidate 
crystal structures. Our numerical scheme elegantly and 
systematically reduces simulation studies of these com- 
plex systems in a two-fold way. Firstly, the problem 
of determining overlaps between irregular particles is re- 
duced tremendously, especially when these particles are 
faceted [4 , by approximating the particle shape with tri- 
angles [24 , since triangle intersections can be straightfor- 
wardly determined. Secondly, predicting candidate crys- 
tal structures is greatly simplified by using the Floppy 
Box Monte Carlo method [25] (FBMC), which efficiently 
analyzes the regular structure on a unit-cell level. We an- 
ticipate our investigation to be the starting point of an 
exploration of (dis) ordered packings of arbitrarily shaped 
particles, which rivals the current research effort on more 
conventional bodies. 

To prove the accuracy of our method, we have repro- 
duced the lower bound ^ib both for the Platonic and for 
the Archimedean Solids, see Ref. [26]. We find excellent 
agreement with Refs. [T5l|l9]-[2T]. The deviation between 
(/)ib and the literature value 0^*^ ranges from 2 • 10 ~^ to 
3 • 10"^. For all results obtained by FBMC (/)ib < (/>rb, be- 
cause compression to the mathematically derived packing 
fraction [151 EI] is virtually impossible. However, FBMC 
performs remarkably well for these convex shapes: the 
simulations we performed typically yielded a very nar- 
row distribution of crystal-structure candidates near the 
closest-packed configuration, the densest of which only 
required minimal additional compression to achieve the 
given value of ^ib- Moreover, for Truncated Tetrahedra 
we have discovered a new crystal structure, a dimer lat- 
tice, with (/)ib = 0.988. . . [26 . This is not only mathe- 



matically interesting, it is also relevant to the study of 
nanoparticle systems, since Truncated Tetrahedral parti- 
cles have recently been synthesized jTI. 

For all prisms and antiprisms with regular n-gonal 
bases (n > 3), for all Johnson Solids, and for all Cata- 
lan Solids, we have verified Ulam's conjecture [27], which 
states that all convex objects pack denser than spheres. 
More importantly, for all of these sets we also find that 
the obtained densest packing for centrally symmetric 
polyhedra is a Bravais lattice packing, in accordance 
with the conjecture of Ref. p^, and that for the con- 
vex, congruent solids without central symmetry it is not 
a Bravais lattice, in accordance with the conjecture of 
Ref. [20 . Here, we present only the data on the Cata- 
lan Solids, which are dual to the Archimedean Solids. In 
Fig. [l] several bounds to the densest packing fraction for 
the Catalan Solids are shown. The outscribed-sphere ap- 
proximation (OSA) and oriented-bounding-box approx- 
imation (OBBA) yield lower bounds to the maximum 
packing fraction, the inscribed-sphere approximation [19] 
(ISA) provides an upper bound. Ref. [26] gives a more 
comprehensive account of the terminology and the em- 
ployed techniques. It also lists the values given in Fig. [T] 
with 6 digit precision. The ISA in combination with the 
FBMC result gives a much narrower region in which the 
densest packing can be found for this group of solids than 
was previously established. The narrowness of the region 
can be attributed to the relatively low asphericity of the 
Catalan Solids. 

Other particle shapes for which we have determined 
the densest-packing crystal structures include the regu- 
lar Pent a- and Heptaprism: ^ib = 0.921 .. . and ^ib = 
0.892 . . . , respectively [26]. The interest in Pentaprisms 
stems not only from the fact that such particles can nowa- 
days be synthesized [28 , but also from the fact that Pen- 
tagons offer the exciting possibility of a structure with 
quasicrystalline order in two dimensions [29 . Hence, 
we consider the prismatic equivalent of Pentagons, Pen- 
taprisms, in our three dimensional (3D) method. Hep- 
taprisms are studied for the purpose of comparison. The 
structures we obtain by FBMC for these Penta- and Hep- 
taprisms are the 3D columnar continuation of the two di- 
mensional (2D) double-lattices found for Pentagons and 
Heptagons, respectively. The prisms in different columns, 
which themselves form a 2D crystal, are out of register. 
The 2D results, to which we compare our 3D crystals, 
were established by mathematical techniques [14^ and 
simulations [18]. Our results indicate that a 3D colum- 
nar continuation of a Pentaprism quasicrystalline con- 
figuration may also exist, since columnar continuations 
appear to be common for non-space-filling prismatic sys- 
tems. For space-filling Tri- and Hexaprisms (and Cubes) 
we observe both columnar continuations and randomly 
stacked layers, i.e., the prisms form crystalline planes, 
but the ordering in different planes is not in register. 

We have also applied our method to study non-convex 
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FIG. 1. Upper and lower bounds to the densest pack- 
ing fraction for the 13 Catalan Solids. The graph shows 
the results for the outscribed-sphere approximation (OSA, 
(/)°b, circles), oriented-bounding-box approximation (OBBA, 
(/)°b5 squares), and FBMC method {(j)ih, connected crosses), 
which give lower bounds to the maximum packing fraction, 
and the inscribed-sphere approximation 19 (ISA, 0°^, dia- 
monds), which gives an upper bound to this packing fraction. 
The value of the densest packing for spheres (/)sph is indicated 
by a red line. A visual representation of the Catalan Solids 
is given below the graph, the abbreviations are as follows: 
Pentagonal Hexecontrahedron (PH), Pentagonal Icositetra- 
hedron (PI), Pentakis Dodecahedron (PD), Disdyakis Tria- 
contrahedron (DT), Deltoidal Hexecontrahedron (DH), Dis- 
dyakis Dodecahedon (DD), Deltoidal Icositetrahedron (Dl), 
Triakis Tetrahedron (TT) , Rhombic Triacontrahedron (RT), 
Triakis Icosahedron (TI), Tetrakis Hexahedron (TH), Small 
Triakis Octahedron (ST), and Rhombic Dodecahedron (RD). 
Note the improvement of FBMC with respect to the OSA 
and OBBA lower bound: the maximum packing fraction can 
henceforth be located in the gray area above the FBMC line. 



(irregular) shapes, which may even contain holes, thereby 
going beyond existing studies. Fig. [2] shows represent a- 




FIG. 2. Four non-convex particle species and their densest-packed crystals. An Echidnahedron (a), a Great Stellated 
Dodecahedron (b), Jessen's Orthogonal Icosahedron (c), and a Stanford Bunny model (d). The crystal structures found by 
FBMC are given in (e) - (h) respectively which show the unit cell and 7 of its periodic images. From left to right the packing 
fractions are (j)\h = 0.294, 0.889, 0.749, and 0.661, for N = 1, 2, 1, and 2 particles in the unit cell, respectively. Note the high 
degree of interpenetration for the Great Stellated Dodecahedron structure, which achieves a surprisingly high packing fraction. 



tions of the shape and predicted crystal structure for 4 
different particles. Additional information on these and 
nine other non-convex bodies can be found in Ref. [26] . 
These candidate crystal structures can be used in theo- 
retical approaches or in simulations with larger system 
sizes to determine their stability using, e.g., free-energy 
calculations. Non-convex particles can pack remarkably 
densely, even when they are not a priori designed to do 
so, as in the case of Escher's Solid. A total of 12 out of 
the 13 shapes achieve a regular packing with ^ib > 0.45, 
and 10 of them achieve (j)\\) > 0.55. The ISA upper bound 
is not very strong in these cases, since we obtain ^°^ = 1 
for all but one of the models, due to their high aspheric- 
ity It is also clear from our result that the packing of 
non-convex (irregular) objects does not possess an intrin- 
sic hierarchy, e.g., in terms of asphericity and symmetry. 

A large fraction of the shapes we have considered is rel- 
atively simple, i.e., the tessellation of the particles only 
involves several dozen triangles. The simulations for such 
particles typically take between 10 minutes and 2 hours 
on a modern 2.0 GHz desktop PC. However, by the gen- 
eral applicability of triangular tessellation, it is also possi- 
ble to study far more complex models. Two such models 
are the Hammerhead Shark, see Ref. [26 , and the Stan- 
ford Bunny, see Fig. [2]i and[2]i, which contain 5,116 and 
3,756 triangles, respectively. Even for these high- triangle 
models the method proves to be relatively quick, the total 
run length of the simulations did not exceed 175 hours. 

The FBMC technique is also perfectly suited to study 
nanoparticle and colloidal systems. As examples we con- 
sider several non-convex shapes of experimentally avail- 
able particles: Caps [5|, Nanostars [3], Tetrapods ^, and 
Octapods [^. Fig. Is] shows representations of these parti- 
cles and some of their crystal structures. Here, we focus 
on the excluded volume interactions due to the complex 



particle shape and neglect any other interactions such 
as. Coulomb, magnetic, and VanderWaals forces. How- 
ever, any additional soft interaction term required to ac- 
curately describe an experimental system can easily be 
included. Even particle deformation is in principle possi- 
ble, but its implementation into the FBMC scheme goes 
beyond the scope of our investigation. 

One physical system we consider in more detail is that 
of the colloidal Cap. The shape of this Cap model, 
see Fig. |3^, is derived from a numerical analysis of the 
collapse of a spherical shell with bending and in-plane 
stretching elasticity submitted to an external isotropic 
pressure [5]. The model we use here is therefore more 
representative of an actual cap than the idealized model 
of a bowl in Ref. [30]. However, the present Cap model 
yields crystal structures similar to the ones described 
in Ref. ^30^, namely: columnar, braided, and inverse 
braided [26]. Fig.pe shows a braided configuration. Sup- 
port is thus provided for the idea that the bowl shape 
captures the essential shape-related physics of these col- 
loidal Cap systems. It should be noted that our simu- 
lations are much more computationally expensive than 
those of Ref. [30], since the Cap tessellation contains 
3,850 triangles in order to obtain a reasonably smooth 
approximation to the actual shape. Typical simulation 
times for Ref. [30] are in the order of 10 minutes, whereas 
ours are around 150 hours. Nevertheless, simulations of 
such complex shapes can be performed with the intent 
to verify the applicability of results obtained for simpler 
models. For faceted nanoparticle and colloid shapes [4], 
as well as smooth objects approximated by a lower num- 
ber of triangles, the simulation times are not prohibitive; 
Nanostars, Octapods, and Tetrapods take the order of 30 
minutes to simulate. 

In conclusion, we have developed a novel method to 
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FIG. 3. Four models of particles which can be synthesized and their densest packing, (a) Two views of a Cap, (b) 
a Nanostar, (c) a Tetrapod, and (d) an Octapod. The structures of the densest-known regular packing of such particles are 
given in (e) - (h) respectively, showing the unit cell and 7 of its periodic images. 



study systems consisting of irregular, non-convex, and 
punctured objects. We have employed this technique to 
determine the densest regular packing and to predict can- 
didate crystal structures in a rigorous way. The complex 
problem of the packing of such shapes has been reduced 
to determining a suitable approximation of a given par- 
ticle in terms of triangles. Using this technique, we have 
predicted candidate crystal structures for non-convex 
particles and we have improved upon the literature values 
for the densest packings of convex solids. In addition, our 
method can easily be extended to study dense amorphous 
(granular) and quasicrystalline packings and systems of 
arbitrarily shaped colloids and nanoparticles in- and out- 
of-equilibrium. Simulations can be performed with or 
without external fields, e.g., gravitational or electric in 
nature, complex inter-particle interaction potentials, or 
even particle shape deformation. Our method thus opens 
the way to a more comprehensive study of the material 
and structure properties than has previously been con- 
sidered feasible. 
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